Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre.

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Microbial eukaryotic predation pressure and biomass at deep-sea hydrothermal vents.

Deep-sea hydrothermal vent geochemistry shapes the foundation of the microbial food web by fueling chemolithoautotrophic microbial activity. Microbial eukaryotes (or protists) play a critical role in hydrothermal vent food webs as consumers and hosts of symbiotic bacteria, and as a nutritional source to higher trophic levels. We measured microbial eukaryotic cell abundance and predation pressure in low-temperature diffuse hydrothermal fluids at the Von Damm and Piccard vent fields along the Mid-Cayman Rise in the Western Caribbean Sea. We present findings from experiments performed under in situ pressure that show cell abundances and grazing rates higher than those done at 1 atmosphere (shipboard ambient pressure); this trend was attributed to the impact of depressurization on cell integrity. A relationship between the protistan grazing rate, prey cell abundance, and temperature of end-member hydrothermal vent fluid was observed at both vent fields, regardless of experimental approach. Our results show substantial protistan biomass at hydrothermally fueled microbial food webs, and when coupled with improved grazing estimates, suggest an important contribution of grazers to the local carbon export and supply of nutrient resources to the deep ocean.

Protistan community composition and metabolism in the North Pacific Subtropical Gyre: Influences of mesoscale eddies and depth.

Marine protists and their metabolic activities are intricately tied to the cycling of nutrients and the flow of energy through microbial food webs. Physiochemical changes in the environment, such as those that result from mesoscale eddies, may impact protistan communities, but the effects that such changes have on protists are poorly known. A metatranscriptomic study was conducted to investigate how eddies affected protists at adjacent cyclonic and anticyclonic eddy sites in the oligotrophic ocean at four depths from 25 to 250 m. Eddy polarity impacted protists at all depths sampled, although the effects of eddy polarity were secondary to the impact of depth across the depth range. Eddy-induced vertical shifts in the water column yielded differences in the cyclonic and anticyclonic eddy protistan communities, and these differences were the most pronounced at and just below the deep chlorophyll maximum. An analysis of transcripts associated with protistan nutritional physiology at 150 m revealed that cyclonic eddies may support a more heterotrophic community, while anticyclonic eddies promote a more phototrophic community. The results of this study indicate that eddies alter the metabolism of protists particularly in the lower euphotic zone and may therefore impact carbon export from the euphotic zone.

Physiology governing diatom vs. dinoflagellate bloom and decline in coastal Santa Monica Bay.

Algal blooms on the Southern California coast are typically dominated by diatom and dinoflagellate taxa, and are governed by their physiological responses to environmental cues; however, we lack a predictive understanding of the environmental controls underlying the establishment and persistence of these distinct bloom events. In this study, we examined gene expression among the numerically dominant diatom and dinoflagellate taxa during spring upwelling bloom events to compare the physiological underpinnings of diatom vs. dinoflagellate bloom dynamics. Diatoms, which bloomed following upwelling events, expressed genes related to dissolved inorganic nitrogen utilization, and genes related to the catabolism of chitin that may have prolonged their bloom duration following nitrogen depletion. Conversely, dinoflagellates bloomed under depleted inorganic nitrogen conditions, exhibited less variation in transcriptional activity, and expressed few genes associated with dissolved inorganic nutrients during their bloom. Dinoflagellate profiles exhibited evidence of proteolysis and heterotrophy that may have enabled them to bloom to high abundances under depleted inorganic nutrients. Taken together, diatom and dinoflagellate transcriptional profiles illustrated guild-specific physiologies that are tuned to respond to and thrive under distinct environmental "windows of opportunity."

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton.

Metagenomics is a powerful method for interpreting the ecological roles and physiological capabilities of mixed microbial communities. Yet, many tools for processing metagenomic data are neither designed to consider eukaryotes nor are they built for an increasing amount of sequence data. EukHeist is an automated pipeline to retrieve eukaryotic and prokaryotic metagenome-assembled genomes (MAGs) from large-scale metagenomic sequence data sets. We developed the EukHeist workflow to specifically process large amounts of both metagenomic and/or metatranscriptomic sequence data in an automated and reproducible fashion. Here, we applied EukHeist to the large-size fraction data (0.8-2,000 µm) from Tara Oceans to recover both eukaryotic and prokaryotic MAGs, which we refer to as TOPAZ (Tara Oceans Particle-Associated MAGs). The TOPAZ MAGs consisted of >900 environmentally relevant eukaryotic MAGs and >4,000 bacterial and archaeal MAGs. The bacterial and archaeal TOPAZ MAGs expand upon the phylogenetic diversity of likely particle- and host-associated taxa. We use these MAGs to demonstrate an approach to infer the putative trophic mode of the recovered eukaryotic MAGs. We also identify ecological cohorts of co-occurring MAGs, which are driven by specific environmental factors and putative host-microbe associations. These data together add to a number of growing resources of environmentally relevant eukaryotic genomic information. Complementary and expanded databases of MAGs, such as those provided through scalable pipelines like EukHeist, stand to advance our understanding of eukaryotic diversity through increased coverage of genomic representatives across the tree of life.IMPORTANCESingle-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

Globally-distributed microbial eukaryotes exhibit endemism at deep-sea hydrothermal vents.

Single-celled microbial eukaryotes inhabit deep-sea hydrothermal vent environments and play critical ecological roles in the vent-associated microbial food web. 18S rRNA amplicon sequencing of diffuse venting fluids from four geographically- and geochemically-distinct hydrothermal vent fields was applied to investigate community diversity patterns among protistan assemblages. The four vent fields include Axial Seamount at the Juan de Fuca Ridge, Sea Cliff and Apollo at the Gorda Ridge, all in the NE Pacific Ocean, and Piccard and Von Damm at the Mid-Cayman Rise in the Caribbean Sea. We describe species diversity patterns with respect to hydrothermal vent field and sample type, identify putative vent endemic microbial eukaryotes, and test how vent fluid geochemistry may influence microbial community diversity. At a semi-global scale, microbial eukaryotic communities at deep-sea vents were composed of similar proportions of dinoflagellates, ciliates, Rhizaria, and stramenopiles. Individual vent fields supported distinct and highly diverse assemblages of protists that included potentially endemic or novel vent-associated strains. These findings represent a census of deep-sea hydrothermal vent protistan communities. Protistan diversity, which is shaped by the hydrothermal vent environment at a local scale, ultimately influences the vent-associated microbial food web and the broader deep-sea carbon cycle.

Methodological 'lenses' influence the characterization of phytoplankton dynamics in a coastal upwelling ecosystem.

New technologies enable the opportunity to improve our monitoring and understanding of marine phytoplankton communities. However, careful consideration for how different methodological approaches, or 'lenses', influence our interpretation of phytoplankton ecology is important, particularly when drawing conclusions about change over time or space. Using both high-throughput 18S rRNA gene sequencing and microscopy, we explored how phytoplankton community structure varied over the course of a year within a nearshore semi-enclosed coastal embayment along the Central Coast of California. The seasonal shift in the relative community dominance (i.e., diatom vs. dinoflagellate dominance) was captured in the microscopy results but not effectively captured in the molecular-based findings. However, the molecular approach explained more of the variability in composition across seasons relative to the microscopy approach. Temporal dynamics of specific bloom-forming taxa also differed between the molecular and microscopy results. Overall, the observed differences between the molecular- and microscopy-derived characterization of phytoplankton dynamics suggest that the approaches are best suited to answer different research questions. Moreover, the approaches complement each other for a more comprehensive perspective of a coastal phytoplankton ecosystem. Therefore, identifying the biases of each approach within natural communities is necessary to effectively and accurately characterize phytoplankton communities.

Daily dynamics of contrasting spring algal blooms in Santa Monica Bay (central Southern California Bight).

Protistan algae (phytoplankton) dominate coastal upwelling ecosystems where they form massive blooms that support the world's most important fisheries and constitute an important sink for atmospheric CO2 . Bloom initiation is well understood, but the biotic and abiotic forces that shape short-term dynamics in community composition are still poorly characterized. Here, high-frequency (daily) changes in relative abundance dynamics of the metabolically active protistan community were followed via expressed 18S V4 rRNA genes (RNA) throughout two algal blooms during the spring of 2018 and 2019 in Santa Monica Bay (central Southern California Bight). A diatom bloom formed after wind-driven, nutrient upwelling events in both years, but different taxa dominated each year. Whereas diatoms bloomed following elevated nutrients and declined after depletion each year, a massive dinoflagellate bloom manifested under relatively low inorganic nitrogen conditions following diatom bloom senescence in 2019 but not 2018. Network analysis revealed associations between diatoms and cercozoan putative parasitic taxa and syndinean parasites during 2019 that may have influenced the demise of the diatoms, and the transition to a dinoflagellate-dominated bloom.

Protistan grazing impacts microbial communities and carbon cycling at deep-sea hydrothermal vents.

Microbial eukaryotes (or protists) in marine ecosystems are a link between primary producers and all higher trophic levels, and the rate at which heterotrophic protistan grazers consume microbial prey is a key mechanism for carbon transport and recycling in microbial food webs. At deep-sea hydrothermal vents, chemosynthetic bacteria and archaea form the base of a food web that functions in the absence of sunlight, but the role of protistan grazers in these highly productive ecosystems is largely unexplored. Here, we pair grazing experiments with a molecular survey to quantify protistan grazing and to characterize the composition of vent-associated protists in low-temperature diffuse venting fluids from Gorda Ridge in the northeast Pacific Ocean. Results reveal protists exert higher predation pressure at vents compared to the surrounding deep seawater environment and may account for consuming 28 to 62% of the daily stock of prokaryotic biomass within discharging hydrothermal vent fluids. The vent-associated protistan community was more species rich relative to the background deep sea, and patterns in the distribution and co-occurrence of vent microbes provide additional insights into potential predator-prey interactions. Ciliates, followed by dinoflagellates, Syndiniales, rhizaria, and stramenopiles, dominated the vent protistan community and included bacterivorous species, species known to host symbionts, and parasites. Our findings provide an estimate of protistan grazing pressure within hydrothermal vent food webs, highlighting the important role that diverse protistan communities play in deep-sea carbon cycling.

Diel transcriptional oscillations of light-sensitive regulatory elements in open-ocean eukaryotic plankton communities.

The 24-h cycle of light and darkness governs daily rhythms of complex behaviors across all domains of life. Intracellular photoreceptors sense specific wavelengths of light that can reset the internal circadian clock and/or elicit distinct phenotypic responses. In the surface ocean, microbial communities additionally modulate nonrhythmic changes in light quality and quantity as they are mixed to different depths. Here, we show that eukaryotic plankton in the North Pacific Subtropical Gyre transcribe genes encoding light-sensitive proteins that may serve as light-activated transcription factors, elicit light-driven electrical/chemical cascades, or initiate secondary messenger-signaling cascades. Overall, the protistan community relies on blue light-sensitive photoreceptors of the cryptochrome/photolyase family, and proteins containing the Light-Oxygen-Voltage (LOV) domain. The greatest diversification occurred within Haptophyta and photosynthetic stramenopiles where the LOV domain was combined with different DNA-binding domains and secondary signal-transduction motifs. Flagellated protists utilize green-light sensory rhodopsins and blue-light helmchromes, potentially underlying phototactic/photophobic and other behaviors toward specific wavelengths of light. Photoreceptors such as phytochromes appear to play minor roles in the North Pacific Subtropical Gyre. Transcript abundance of environmental light-sensitive protein-encoding genes that display diel patterns are found to primarily peak at dawn. The exceptions are the LOV-domain transcription factors with peaks in transcript abundances at different times and putative phototaxis photoreceptors transcribed throughout the day. Together, these data illustrate the diversity of light-sensitive proteins that may allow disparate groups of protists to respond to light and potentially synchronize patterns of growth, division, and mortality within the dynamic ocean environment.

A Primer for Microbiome Time-Series Analysis.

Time-series can provide critical insights into the structure and function of microbial communities. The analysis of temporal data warrants statistical considerations, distinct from comparative microbiome studies, to address ecological questions. This primer identifies unique challenges and approaches for analyzing microbiome time-series. In doing so, we focus on (1) identifying compositionally similar samples, (2) inferring putative interactions among populations, and (3) detecting periodic signals. We connect theory, code and data via a series of hands-on modules with a motivating biological question centered on marine microbial ecology. The topics of the modules include characterizing shifts in community structure and activity, identifying expression levels with a diel periodic signal, and identifying putative interactions within a complex community. Modules are presented as self-contained, open-access, interactive tutorials in R and Matlab. Throughout, we highlight statistical considerations for dealing with autocorrelated and compositional data, with an eye to improving the robustness of inferences from microbiome time-series. In doing so, we hope that this primer helps to broaden the use of time-series analytic methods within the microbial ecology research community.

Biological composition and microbial dynamics of sinking particulate organic matter at abyssal depths in the oligotrophic open ocean.

Sinking particles are a critical conduit for the export of organic material from surface waters to the deep ocean. Despite their importance in oceanic carbon cycling and export, little is known about the biotic composition, origins, and variability of sinking particles reaching abyssal depths. Here, we analyzed particle-associated nucleic acids captured and preserved in sediment traps at 4,000-m depth in the North Pacific Subtropical Gyre. Over the 9-month time-series, Bacteria dominated both the rRNA-gene and rRNA pools, followed by eukaryotes (protists and animals) and trace amounts of Archaea. Deep-sea piezophile-like Gammaproteobacteria, along with Epsilonproteobacteria, comprised >80% of the bacterial inventory. Protists (mostly Rhizaria, Syndinales, and ciliates) and metazoa (predominantly pelagic mollusks and cnidarians) were the most common sinking particle-associated eukaryotes. Some near-surface water-derived eukaryotes, especially Foraminifera, Radiolaria, and pteropods, varied greatly in their abundance patterns, presumably due to sporadic export events. The dominance of piezophile-like Gammaproteobacteria and Epsilonproteobacteria, along with the prevalence of their nitrogen cycling-associated gene transcripts, suggested a central role for these bacteria in the mineralization and biogeochemical transformation of sinking particulate organic matter in the deep ocean. Our data also reflected several different modes of particle export dynamics, including summer export, more stochastic inputs from the upper water column by protists and pteropods, and contributions from sinking mid- and deep-water organisms. In total, our observations revealed the variable and heterogeneous biological origins and microbial activities of sinking particles that connect their downward transport, transformation, and degradation to deep-sea biogeochemical processes.

SSU-rRNA Gene Sequencing Survey of Benthic Microbial Eukaryotes from Guaymas Basin Hydrothermal Vent.

Microbial eukaryotes have important roles in marine food webs, but their diversity and activities in hydrothermal vent ecosystems are poorly characterized. In this study, we analyzed microbial eukaryotic communities associated with bacterial (Beggiatoa) mats in the 2,000 m deep-sea Guaymas Basin hydrothermal vent system using 18S rRNA gene high-throughput sequencing of the V4 region. We detected 6,954 distinct Operational Taxonomic Units (OTUs) across various mat systems. Of the sequences that aligned with known protistan phylotypes, most were affiliated with alveolates (especially dinoflagellates and ciliates) and cercozoans. OTU richness and community structure differed among sediment habitats (e.g. different mat types and cold sediments away from mats). Additionally, full-length 18S rRNA genes amplified and cloned from single cells revealed the identities of some of the most commonly encountered, active ciliates in this hydrothermal vent ecosystem. Observations and experiments were also conducted to demonstrate that ciliates were trophically active and ingesting fluorescent bacteria or Beggiatoa trichomes. Our work suggests that the active and diverse protistan community at the Guaymas Basin hydrothermal vent ecosystem likely consumes substantial amounts of bacterial biomass, and that the different habitats, often defined by distances of just a few 10s of cm, select for particular assemblages and levels of diversity.

Are We Overestimating Protistan Diversity in Nature?

Documenting the immense diversity of single-celled, eukaryotic organisms (protists) has been a formidable challenge for ecologists. These species were originally defined by morphological criteria, but shortcomings of the morphospecies concept, and a bewildering array of sizes and cellular attributes, has made constructing a taxonomy that is useful for ecologists nearly impossible. Consequently, physiological and genetic information has been integrated to address these shortcomings, and to develop the framework of a unifying taxonomy. DNA sequence information, in particular, has revolutionized studies of protistan diversity. However, the exponential increase in sequence-based protistan species richness published from field surveys in recent years raises the question of whether we have moved beyond characterizing species-level diversity and begun to reveal intraspecies diversity. The answer to that question appears to be 'yes', at least for some protistan lineages. The need to document such microdiversity may be justified, but it is important for protistologists to recognize and acknowledge that possibility, and its consequences.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone.

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

Single-cell transcriptomics of small microbial eukaryotes: limitations and potential.

Single-cell transcriptomics is an emerging research tool that has huge untapped potential in the study of microbial eukaryotes. Its application has been tested in microbial eukaryotes 50 µm or larger, and it generated transcriptomes similar to those obtained from culture-based RNA-seq. However, microbial eukaryotes have a wide range of sizes and can be as small as 1 µm. Single-cell RNA-seq was tested in two smaller protists (8 and 15 µm). Transcript recovery rate was much lower and randomness in observed gene expression levels was much higher in single-cell transcriptomes than those derived from bulk cultures of cells. We found that the reason of such observation is that the smaller organisms had much lower mRNA copy numbers. We discuss the application of single-cell RNA-seq in studying smaller microbial eukaryotes in the context of these limitations.

Protistan diversity and activity inferred from RNA and DNA at a coastal ocean site in the eastern North Pacific.

Microbial eukaryotes fulfill key ecological positions in marine food webs. Molecular approaches that connect protistan diversity and biogeography to their diverse metabolisms will greatly improve our understanding of marine ecosystem function. The majority of molecular-based studies to date use 18S rRNA gene sequencing to characterize natural microbial assemblages, but this approach does not necessarily discriminate between active and non-active cells. We incorporated RNA sequencing into standard 18S rRNA gene sequence surveys with the purpose of assessing those members of the protistan community contributing to biogeochemical cycling (active organisms), using the ratio of cDNA (reverse transcribed from total RNA) to 18S rRNA gene sequences within major protistan taxonomic groups. Trophically important phytoplankton, such as diatoms and chlorophytes exhibited seasonal trends in relative activity. Additionally, both radiolaria and ciliates displayed previously unreported high relative activities below the euphotic zone. This study sheds new light on the relative metabolic activity of specific protistan groups and how microbial communities respond to changing environmental conditions.

Estimating Protistan Diversity Using High-Throughput Sequencing.

Sequencing hypervariable regions from the 18S rRNA gene is commonly employed to characterize protistan biodiversity, yet there are concerns that short reads do not provide the same taxonomic resolution as full-length sequences. A total of 7,432 full-length sequences were used to perform an in silico analysis of how sequences of various lengths and target regions impact downstream ecological interpretations. Sequences that were longer than 400 nucleotides and included the V4 hypervariable region generated results similar to those derived from full-length 18S rRNA gene sequences. Present high-throughput sequencing capabilities are approaching protistan diversity estimation comparable to whole gene sequences.

Investigating microbial eukaryotic diversity from a global census: insights from a comparison of pyrotag and full-length sequences of 18S rRNA genes.

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.